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ABSTRACT 

We present a ~6.5' x8' Expanded Very Large Array (EVLA) mosaic observations of the NH3 (1,1) emission 
in the Barnard 5 region in Perseus, with an angular resolution of 6". This map covers the coherent region, 
where the dense gas presents subsonic non-thermal motions (as seen from single dish observations with the 
Green Bank Telescope, GBT). The combined EVLA and GBT observations reveal, for the first time, a striking 
filamentary structure (20" wide or 5,000 AU at the distance of Perseus) in this low-mass star forming region. 
The integrated intensity profile of this structure is consistent with models of an isothermal filament in hydro- 
static equilibrium. The observed separation between the B5-IRS1 young stellar object (YSO), in the central 
region of the core, and the northern starless condensation matches the Jeans length of the dense gas. This sug- 
gests that the dense gas in the coherent region is fragmenting. The region observed displays a narrow velocity 
dispersion, where most of the gas shows evidence for subsonic turbulence, and where little spatial variations 
are present. It is only close to the YSO where an increase in the velocity dispersion is found, but still displaying 
subsonic non-thermal motions. 

Subject headings: ISM: clouds — stars: formation — ISM: molecules — ISM: individual (Perseus Molecular 
Complex, B5) 



1. INTRODUCTION 

Molecular clouds (MCs) observed using low-density trac- 
ers display velocity dispersions much larger than the thermal 
value (e.g., |Zuckerman & Evans 1974; Larson 1981; Myers 
1983). These supersonic motions are usually attributed to 
"turbulence," and a variety of recent numerical models of tur- 
bulence can reproduce qualitatively realistic clouds. Measure- 
ments of MCs' energy budgets show that turbulence must be 
dissipated in order for dense cores (where individual or small 
groups of stars form) to collapse and form stars. 

Dense cores have been studied using NH 3 (1,1) maps, 
which traces material with densities of a few 10 4 cm -3 . It 
has been found that dense cores show an almost constant 
level o f non-thermal motions , within a certain "coherence" 
radius (Good man et al.|[T998| ). The term "coherent core" is 
used to describe the dense gas where non-thermal motions 
are roughly constant, and typically small er than the therm al 
motions, independent of scale (see also Caselr Tet" al.|2002| ). 

Recent observations carried out with the 100-m Robert F. 
Byrd Green Bank Telescope (GBT) have allowed us to cre- 
ate large scale maps of NH 3 (1,1) towards four star-forming 
regions in the Per seus Molecular Cloud: B5, IC348-SW, 
L1448, a ndL1451 ([Pineda et al.|2QT0l[2QTT| ). One of the main 
results of Pineda et al. ( 2010]) is the clear observation (for the 
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first time) of the sharp transition (at 30" resolution) between 
the coherent section of the B5 dense core and the mor e tur- 
bulent dense gas outside it (see also Pineda et al.||20lT| for 
results on the other regions surveyed). However, the angular 
resolution of these observations did not allow us to study in 
great detail the spatial variations of the velocity dispersion or 
column density. 

In this letter, we present new NH3 (1,1) observations of 
B5 obtained with the Expanded Very Large Array (EVLA) 
throu gh the Open Sha red Risk Observing (OSRO) program, 
see |Perley et al.| ( |20TT] ) for a description of the EVLA project. 
These data are combined with previous observations from the 
GBT to study the dense gas traced by NH 3 at high angular and 
spectral resolution. Here, we present results from the com- 
bined EVLA and GBT observations that enable us to analyze 
the kinematic properties of the dense gas traced by NH 3 (1,1) 
and also to study the radial profile of new filamentary struc- 
ture found within this core. 

2. DATA 

Observations were carried out with the EVLA of the Na- 
tional Radio Astronomy Observatory on May 16th, 2010 
(project 10A-181). We used the high-frequency K-band 
receiver and configured the WIDAR correlator to observe 
a 1 MHz window around the NH3 (1,1) rest frequency 
(23.6944950787 GHz). The correlator generated 256 chan- 
nels across the window, giving a 3.90625 kHz channel sepa- 
ration, equivalent to 0.049 krns" 1 at the observed frequency. 
This configuration covers the main hyperfine component and 
also one of the inner pair of satellite lines for NH3 (1,1). At 
this frequency the primary beam of the array is about 1.9'. 
The array was in the compact (D) configuration, which covers 
baselines from 35 meters up to 1 km, and included 26 anten- 
nas. The observations covered the entire region of interest 
using 27 pointings. The observations were carried out under 
fair weather (T22GHZ = 0.0689 nepers) for 7 hours. 

The bandpass and absolute flux calibrator was the quasar 
3C48, while the phase and amplitude calibrator, quasar 
J0336+3218, was observed every 15 minutes. The X-band 
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reference pointing checks were performed every 60 minutes 
using the quasar J0336+3218. 

The raw-data were reduced using CASA image processing 
software. The signal from each baseline was inspected, and 
baselines showing spurious data were removed prior to imag- 
ing. Each channel was cleaned separately according to the 
spatial distribution of the emission. The clean map was im- 
aged using a 6" circular beam and corrected by the primary 
beam. The images were created using multi- scale clean (with 
scales of 4 and 12 arcsec and smallscalebias=0.2) with a ro- 
bust parameter of 0.5. 

Imaging of the EVLA data only produced a resulting data 
cube where strong negative bowls were found. These nega- 
tive bowls are generated by the large scale NH 3 (1,1) emission 
which arises from angular scales not sampled by the interfer- 
ometer. We in cluded the NH3 (1 1) single dish data obtained 
with the GBT ( Pineda et al. 2010) as a prior (or model) during 
the imaging process to recover the large scale emission. The 
final noise achieved is 14mJybeam _1 per channel. 

The integrated intensity map (see Figure [T]) shows a few re- 
gions with negative emission, which wouldsuggest the pres- 
ence of imaging artifacts. However, careful inspection of the 
data confirms that the spectra towards those positions do not 
show any clear negative bowl associated with imaging arti- 
facts, in fact, these regions are about the -5-a level still con- 
sistent with the map noise. 



3. RESULTS 

Left panel of Figure [T] shows the integrated intensity map 
obtained using the GBT at 30" resolution. Gray contours in 
Figure [T] show the extent of NH3 (1,1) emission. The orange 
contours show the regions in the GBT data within which the 
non-thermal velocity dispersion is sub-sonic. The blue con- 
tour shows the region observed with the EVLA and presented 
in the right panel, and it covers the entire central region where 
the sub-sonic non-thermal velocity dispersions are observed. 
The resulting NH 3 (1,1) integrated intensity map for B5 is 
shown in right panel of Figure [Tl and it covers a region of size 
6.5 / x8 / . 

The new high resolution integrated intensity map (right 
panel of Figure [lj reveals that within the region of subsonic 
non-thermal motions found in the single dish data (orange 
contour) filamentary structures appear. These filaments are 
narrow, with width s of ~ 20" or 5,000 AU at the distance of 
Hirota et al.|2008b. 



Perseus (250 pc; 

An important feature of molecular line observations is the 
ability to probe the kinematics of the gas. Here we fit simulta- 
neously all hyperfine components of the N H3 (1,1) line using 
a forward mode l previously present ed by |Rosolowsky et al. 
(2008), see also Pineda ^taLl ( |20 1 0| ) . This method describes 
the emission at every position with a centroid velocity (vlsr), 
velocity dispersion (cr v ), kinetic temperature (7^), excitation 
temperature (T ex ) and opacity (rn), while also including the 
response of the frequency channel using a sine profile. Since 
the kinetic temperature is only used to predict the NH 3 (2,2) 
line (not observed due to the constraints in the OSRO pro- 
gram), it does not have any effect on the remaining parame- 
ters. We use a fixed value of the kinetic temperature of 10 K 
for the en tire region, which is consistent with the results ob- 
tained by Pin eda et al.| (|20TTJ using the single dish data. If 
the resulting fit does not provide an accurate velocity disper- 
sion, (T v < 0.05 km s -1 (velocity dispersion narrower than the 
expected thermal value for gas at 6 K) or o Gv > 0.2 a v (signal- 



to-noise for the velocity dispersion less than 5), the fit is re- 
peated but with a fixed value of 5 K for excitation temperature. 
The centroid velocity and velocity dispersion maps are shown 
in Figure [2] 

The centroid velocity map shows little variation (< 
0.6 km s -1 ) across the entire region, see Figure^. The region 
close to the YSO displays a velocity gradient at a position 
angle that is close to the outflow (shown by the arrows) but 
in opposite direction. The velocity dispersion map, see Fig- 
ure shows vast regions where small and uniform velocity 
dispersion are found. It is only towards positions close to the 
YSO where slightly broader lines are found. 

Figure [3] presents the distribution of the derived velocit 
dispersion towards B5. Two histograms are shown Figure r 
depending on the proximity to the YSO in B5: (a) positions 
close to the YSO (separated by 2 beams or less, <12") in red, 
and (b) all other pixels in black. For points far from the YSO 
the velocity dispersions are small and the velocity dispersion 
distribution (black histogram in Figure [3]) is narrow, almost 
every pixel at a distance larger than 2 beams (12") from the 
central YSO displays sub-sonic non-thermal motions. The 
velocity dispersion distribution of pixels far from the YSO, 
black histogram in Figure [3] peaks at a value lower than what 
is expected if the non-thermal component, a^r, is equal to 
half the thermal velocity dispersion, 0.5c s ave . The velocity 
dispersion of points close to the YSO, red histogram in Fig- 
urel3| presents lines broader than the rest of the core, but they 
stilfdisplay velocity dispersions with a subsonic non-thermal 
component. This increase in the velocity dispersion might be 
the effect of the radiation from the embedded YSO or due to 
the interaction between the outflow or stellar wind and the 
dense gas. 

Herschel observations o f the IC 5146 star forming region 
( [Arzoumanian et aL 201 1 ) revealed filamentary structure seen 
in the column d ensity maps (see also Andre et al.|2010b |Ar-| 
zoumanian et a l. (2011) fitted the column density profile of 
filaments with a cylindrical filament model, 
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where p c is the filament's central density, r is the cylindri- 
cal radius, p is the power-law density exponent at large radii, 
Rfi at is the radius of the density profile inner flat section, and 
A p is a finite constant factor dependent on p and the filament 
inclination angle. They find filaments which are wel l fit with a 
densit y exponent in the range p = 1 .5-2.5 (see also |Lada et al.| 
1999] ). These exponents do not agree with the predicted val- 
ues of an isothermal filament in hydros tatic equilibri um where 
a steeper exponent, p = 4, is expected ( Ostrik er|1964| ). 

Here we focus our attention towards the filament shown in 
Figure [Tlby the yellow box. Since this filament is almost per- 
fectly aligned in the North-South direction we perform a se- 
ries of horizontal cuts, and define the radius as the distance 
from the peak at a given cut. The average velocity disper- 
sion and integrated intensity emission profiles along the fil- 
ament are shown in panels (a) and (b) of Figure |4| respec- 
tively, where the spread in the distribution is shown by the 
yellow area. Figure [4^ shows that the velocity dispersion 
does not change across the filament, and it is consistent with 
subsonic non-thermal velocity dispersions (delimited by the 
dotted line). Figure [4J3 shows the integrated intensity profile, 
which is much wider than the beam of the combined EVLA 
and GBT data (shown by the blue dash line). Since we do 
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Figure 1. Left panel: Integrated intensity map of B5 in NH3 (1,1) obtained with GBT. Gray contours show the 0.15 and 0.3 Kkms -1 level in NH3 (1,1) integrated 
intensity. The orange contours show the region in the GBT data where the non-thermal velocity dispersion is sub-sonic. The young star, B5-IRS1, is shown by 
the star in both panels. The outflow direction is shown by the arrows. The blue contour shows the area observed with the EVLA, and the red box shows the 
area shown in the right panel. Right panel: Integrated intensity map of B5 in NH3 (1,1) obtained combining the EVLA and GBT data. Black contour shows the 
50mJybeam _1 kms -1 level in NH3 (1,1) integrated intensity. The yellow box shows the region used in Figure]?] The northern starless condensation is shown by 
the dashed circle. 
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Figure 2. (a) Centroid velocity map of B5 obtained by fitting all observed NH3 hyperfine components simultaneously. The position of the YSO is shown by the 
star, and the orientation of the outflow is shown by the arrows. The estimated Jeans length, Ay, of the dense gas is shown by the vertical scale bar, see Section[4] 
The EVLA beam size is shown at the bottom left. The orange contour shows the region with subsonic non-thermal motions identified from the GBT data. Gray 
contour shows the 50mJybeam _1 kms -1 level in NH3 (1,1) integrated intensity, same as right panel of Figure^ The Black box shows the region used in Figure]?] 
(b) Velocity dispersion map derived from fitting all NH3 hyperfine components simultaneously. 
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Figure 3. Velocity dispersion distribution for positions at a distance to B5- 
IRS1 larger (smaller) than 12" is shown in black (red). Blue vertical lines 
show the expected velocity dispersion for two values of the velocity disper- 
sion non-thermal component (&nt)'- 0.5c s ,ave, and c s ,ave, where c s ,ave is the 
sound speed of the average particle (/i = 2.33) assuming Tki n =10 K. The ve- 
locity dispersions, cr v , obtained far enough from the YSO (black histogram) 
are consistently below the values expected for supersonic non-thermal mo- 
tions. Positions close to the YSO (red histogram) display larger velocity dis- 
persions than the rest of the region, but they are still consistent with subsonic 
non-thermal motions. 

not observe the NH 3 (2,2) with the EVLA we cannot provide 
a column density profile, but we use NH3 (1,1) integrated in- 
tensity as a proxy for the total column density ( [Rosolowsky 



et al.[2008[|Friesen et al.|20Q9l poster et al.|2009[jPineda et al.| 
The integrated intensity profile is fitted with a model 



follows equation [T] and the best fit results for p = 2 and 
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p = 4 are shown in Figure [4] The figure show s that this fila- 
ment in B5, in opposition to what is found by Arzoumanian 
|et al.| ( [20TTT ) on larger filaments, is better fit by an isothermal 
hydrostatic filament model (see also Hacar et al., submitted, 
or Bourke et al., submitted, for other examples of objects with 
similar column density profiles). 

4. DISCUSSION AND CONCLUSION 

The observations presented here show that subsonic non- 
thermal velocity dispersions display little variations across 
the region of coherence, and even across a filament. Also, 
the coherent region is far from presenting an uniform column 
density (as traced by the NH3 (1,1) integrated intensity), as 
filamentary substructures are revealed. But the presence of 
substructures within the region of coherence should not be 
surprising, the lack of turbulent support might allow exter- 
nal forces or Jeans-like instabilities to generate over-densities 
without much difficulties. 

Filaments have received special attention recently. They 
appear as a natural ou tcome from turbulen t simulations of 
molecular clouds (e.g., |Klessen et al.||2004]), and a lso from 
the fragmentation of a modulated layer (Myers 2009). Recent 
Herschel observations have shown that filaments are c om- 
monly found in star-for ming regions (Andre et al. 2010|[Ar j 
zoumania n et al.||2011| ). However, these filaments present 
different properties than those presented here : (a) they are 
bigger structures, with a characteristic width of ~0.1 pc; and 
(b) they are not well fitted by an hydrostatic and isothermal 
cylindrical model. The filamentary structures revealed by 
our observations appear to have different properties than fila- 
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Figure 4. (a) Average velocity dispersion profile perpendicular to the fil- 
ament shown by the yellow box in Figure [ll The dispersion of the radial 
profile along the filament is shown by the yellow area. Similar to Figure [3] 
blue horizontal lines show the expected velocity dispersion for two values 
a^T'- 0.5c s ,ave, and c s ,ave- (b) Average integrated emission profile perpendic- 
ular to the filament, with the dispersion shown in yellow. The beam response 
is shown by the blue line, while two filament models are shown by the red 
curves. 



ments previously studied. The filaments in B5 are embedded 
in a dense region with subsonic turbulence, and the filament 
shown in Figure R] is better fitted with the hydrostatic isother- 
mal filament moael. The differences between the filaments 
identified using Herschel and EVLA are due to physical dif- 
ferences between the structures. The Herschel identified fil- 
aments are found at low column densities and have a char- 
acteristic size of ~ 0.1 pc, which is similar to the size of the 
coherent structures found in sing le dish observations dPineda| 
|et al.|2010[|Goodman et al.| 1998 ), but much smaller than the 
largest scales detected (~ 1 deg). While the filaments revealed 
by the EVLA observations are substructures within the coher- 
ent region which are also at higher column densities than the 
structures found in Herschel. Also, since the filament width 
(see Figure HJ)) is smaller than the EVLA primary beam and 
the single dish data is used to recover the extended emission 
then the radial profile is little affected by the spatial filtering 
of interferometric observations. 

Fragmentation in the coherent region should occur at the 
Jeans length, 



A., = 44,300 



10K 



1/2 



(n> 



10 4 cnr 



-1/2 



AU, (2) 



where T is the temperature of the gas, (n) is the average den- 
sity, and fx is the mean molecular weight. For the average 
values of density and temperature in B5 from the dense gas 
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(T = 10 K and (n) = 10 4 cnT 3 , [Pineda et al.|201l} the Jeans 



length correspond to the separation between the B5-IRS1 
YSO and the starless condensation in the northern part of the 
region observed. Fragmentation has been argued to e xplain 
the stars separation in young star-form ing regions (e.g., Hart- 
mann 2002; Teixeira et al. 2006, 2007 \ and the mass function 
of cores or condensation s (e.g. , |Testi & Sargent|19 98; Motte 
|et al.||1998t | Andre et al.||2007| ). Here we show evidence for 
fragmentation occurring within the coherent region, which is 
not necessarily linked to the filament formation process. 

These results support the idea that a coherent region needs 
to be created to form a low-mass star. It is in this coherent re- 
gion where the dense gas that has lost its turbulent support ac- 
cumulates, and then it can easily fragment, create a filament, 
and/or undergo gravitational collapse to form a star. We ex- 
pect that future EVLA and ALMA observations of star form- 
ing regions will allow us to produce temperature and velocity 
dispersion maps which will allow us to perform a much better 
job on characterising the physical properties of filaments and 
of these coherent regions. 
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